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Abstract: We study the single production of charged Higgs bosons in e + e colli- 
sions, chiefly in the Minimal Supersymmetric Standard Model. Our analysis com- 
^ ■ plements foregoing studies of the pair production channel especially in regions where 

the kinematic constraint suppresses pair production. We present cross sections rele- 
vant to experiments at the next generation linear colliders and some brief discussions 
of their phenomenology. Our analysis shows that the single production is a useful 
alternative channel for studying the phenomenology of charged Higgs bosons, and 
there are regions of parameter space where it is accessible beyond the kinematic limit 
for pair production. 
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1. Introduction 

Charged Higgs bosons H ± are a cornerstone of beyond-the-Standard-Model phe- 
nomenology. They prediction of the supersymmetric extensions of the 
Standard Model (SM) from purely theoretical requirements of consistency. Their 
discovery and the confirmation of their properties will be a significant step towards 
a full understanding of electroweak symmetry breaking. 

Some substantial effort has thus been channelled into the evaluation of their 
phenomenology at future colliders, and over the past few years the situation with 
respect to their discovery potential has become more clear. 

At hadron colliders [Q, especially with the luminosity of the LHC, charged Higgs 
bosons below the top quark mass are expected to be produced abundantly in the 
decay of top quarks. When their mass is near or greater than the top quark decay 
threshold, their discovery prospects are further hindered by the falling structure 
function but they will be produced mainly in the 'tfr-fusion' process, namely the 
parton level process gb — > tH~ and the charge conjugate. There have recently been 
attempts [0] to connect the two channels together by looking at a generic subprocess 
gg — > tbH~ which includes a component of the latter process that is leading order at 
large transverse momentum of the 'spectator' bottom quark. In both cases the most 
promising decay mode is presumably the H — > t~v mode ||. It is claimed that 
relatively heavy charged Higgs bosons, of mass up to 1 TeV, can be probed using 
this mode, especially at large tan /3, although a full experimental simulation similar 
to that presented in [|J will be essential to test this claim and to draw the discovery 
contour. Here tan/3 is as usual defined as the ratio of the vacuum expectation 
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values of the two Higgs doublets. Several other production channels have also been 
considered J5|. 

At future electron-positron linear colliders (LCs), the dominant production pro- 
cess, if kinematically allowed, is the pair production process 0: 

e~e + ^H + H-. (1.1) 

The production rate depends only on the charged Higgs mass and the process provides 
a hallmark channel through which we can study H ± phenomenology J/J. 

When the charged Higgs mass is near or above half the centre-of-mass energy the 
phenomenology is much less well understood. There has lately been some interest in 
the associated production mode with W ± [§, || p~0|| : 



e e 



H ± W T , (1.2) 



although the extent to which this channel could contribute towards the study of 
charged Higgs phenomenology is not clear. In models with a multi-Higgs-doublet 
structure, the process is loop-induced in the massless electron limit, and therefore 
the cross section is suppressed. We also note that no analysis of the SM background 
is available so far. 

Other channels that have been studied in the literature are the heavy quark 



associated production mode |TT 



e-e + -> bbW ± H T , (1.3) 
and the tv t associated production mode 



e-e + -> H + r-u T , H-r + u T . (1.4) 

The purpose of this paper is to study the above processes in greater detail and to 
collect them together with many other single production channels of charged Higgs 
bosons in order to complement the pair production channel both above and below its 
kinematic threshold. Although our method is applicable generally in principle, we 
have adopted the Minimal Supersymmetric Standard Model (MSSM) for calculating 
Higgs masses and mixings. The MSSM being a model with a decoupling structure 
in the Higgs sector, the resulting cross sections are small compared to cases where, 
for example, extra resonances are available. However, it is beyond our intention to 
study models other than the MSSM in the present paper. 

2. Production processes 

We consider the following fourteen processes: 

e~e + -> t-v t H + ,t + v t H- (2.1) 
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ing to the above are shown in figure [I]. 





For all tree-level processes, the matrix elements were calculated both by hand, 
for select set of diagrams in order to simplify the calculation, using the usual trace 



method, and by using the helicity amplitude formalism of |13[ . The two sets of results 
agree up to the numerical precision employed. We present the matrix element squared 
for the dominant tree-level signal processes (|2.1|) and ( |2.2| ) in the appendix. 

All processes were calculated at leading order only. For the 2HDM parameters, 
we adopted the MSSM throughout. For the SM parameters we adopted the following: 
m b = 4.25 GeV, m t = 175 GeV, m e = 0.511 MeV, m T = 1.78 GeV, m v = 0, M w = 
80.23 GeV, T w = 2.08 GeV, M z = 91.19 GeV, T z = 2.50 GeV, sin 2 6 W = 0.232. The 
top quark width T t was evaluated at leading order for each value of M H ± and tan (3. 
Neutral and charged Higgs masses were calculated for given values of M40 and tan /3 



using the HDECAY package |L4|, with the SUSY masses, the trilinear couplings and 
the Higgsino mass parameter \x being set to 1 TeV. The charged Higgs boson width 
T H ± was evaluated at leading order. 



The loop- induced processes ( |2~3| ) and ( |2.4| ) were adapted from ||, [L5|]. In the one- 



loop analyses, we assumed that the superpartners are sufficiently heavy to decouple, 
so that only the heavy quark loops and Higgs-gauge loops are included. We intro- 
duced counter terms from WH and wH mixings (w represents the Nambu-Goldstone 
boson) |l6j and used the renormalised coupling to account for the mixing. Details 
of the calculation are shown in ||. One of the renormalisation conditions is that the 
WH mixing is zero for onshell H ± , such that the third diagram of figure |I](b) is ef- 
fectively zero. The cross sections are qualitatively consistent with other calculations 
[§, |K|, although we note that this renormalisation is not performed in |J. 

The HWZ and HW^ vertices, which enter into the process ([2.4D, have been 
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calculated in [15] and [17|] for the HWZ vertex and in (18| for both vertices. In 
our analysis, we improved on their calculations by renormalising the WH and wH 
mixings according to the calculation in and [[H]]. Again, the superpartners were 



assumed to be heavy. In addition, for process Q2.4|), we dropped the bosonic loop 
contributions in order to save time during numerical computation. This does not 
affect the results, as the contributions from Higgs-gauge loop diagrams are small 
in parameter regions where these vertices are substantial. We have explicitly veri- 
fied this statement numerically. The total cross sections were evaluated using these 
vertices with the helicity amplitude method. 

In the results which we present, the charge conjugate subprocesses are not in- 
cluded. Our results correspond to the production of either an H + or an H~ (except 
for the pair production process). This is in order to avoid double counting final 
states. For instance, if we consider process ( |2.1| ), the final state can be t~v t r + u T . 
Far from the kinematic threshold for pair production the 'total' single charged Higgs 
production for a final state XH ± may be given by: 

a{XH ± ) = a(XH + ) + a{XH~) - BR^ X)a(H + H-). (2.15) 

Near the kinematic threshold we can not treat the processes consistently using this 
formula without specifying the particular final state, and this would limit the gener- 
ality of our approach. The total cross section is twice that presented here in regions 
where pair production is forbidden, and nearly equal to that presented here in the 
limit where the branching ratio BR{H ± — > X) tends to one. 



3. Production cross sections 

We present the cross sections as functions of the charged Higgs boson mass M#± 
at collider energies of 500 GeV and 1000 GeV and 4 different values of tan/3, 1.5, 
7, 30 and 40. These are shown in figures [| to || If we assume, for instance, an 
integrated luminosity of 500 fb _1 Jl9|], 10~ 5 pb corresponds to 5 events per year 
before acceptance cuts and background reduction. We do not discuss the background 
reduction procedure in detail in this study, and 10~ 5 pb is taken naively as the 
threshold of the 'relevance' of the process to the study of charged Higgs production at 
LC. We emphasise that this is not intended in any way as a threshold of detectability, 
or even visibility, as the evaluation of such thresholds would require jet simulations 
and machine-dependent considerations which are clearly beyond the scope of this 
current study. 

For the sake of comparison we also present the cross sections of the on-shell pair 
production mode ( |TTT| ) in figure []. 

In the allowed kinematic ranges, the tv t and tb associated production processes 
( |2.1| ) and ( |2.2j ), shown in figure |2|, are dominated by pair production. When M#± is 
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small, process ( |2.2|) also has a large contribution from top pair production followed 
by the decay of one of the top quarks into bH ± , which explains the rise of the 
cross section below 175 GeV. Beyond the kinematic limit for pair production, which 
occurs at M H ± ~ y/s/2, the cross sections still exceed 10 -5 pb for some values of 
tan (3. Process ( |2.1| ) is enhanced for large tan (3 whereas process ( |2.2| ) is enhanced for 
both large and small values of tan j3 and the minimum is at tan f3 = \Jmt/ mb ~ 7. 

Figure || shows the rate for the loop-induced W ± H T associated production pro- 
cess (2J). When tan/5 is small, the top Yukawa coupling is enhanced and the cross 
section is large. Hence this channel complements process (|2.1|) . We note that the 
tan (3 dependence of the signal rate for this process is ~ m\ cot 2 (3 at small tan (3 
and ~ ml tan 2 (3 at very large tan (3. Hence the bottom Yukawa contribution is 
suppressed. The cross section remains large beyond the pair production kinematic 
limit. The peak in the cross section near 200 GeV corresponds to the threshold 
m t + rrib ~ M#± , after which the cross section falls slowly up to the kinematic limit 
at M H ± ~ y/s — Mw± ■ 

We note that at tan/5 ~ 7, both the t~v t H + mode and the W ± H Zf mode have 
cross sections near 10~ 5 pb at \fs = 500 GeV and M H ± > 250 GeV. Thus there 
are always regions in M#± where charged Higgs bosons can be produced at a 500 
GeV machine even when the pair production channel is unavailable. However, we 
should bear in mind that the detectability of such processes needs more realistic sim- 
ulations in order to quantify the effect of the decay, the fragmentation and detector 
resolutions. 

In figure |] we show the rate for the loop induced vector boson fusion process 
( [2.41 ). The overall rate is small at both 500 GeV and 1 TeV, but there is an interesting 
tan (3 dependence of the cross section where, near the tb — > H + kinematic threshold, 
the rates are enhanced for large tan/3. As with the other vector fusion induced 
processes which we present later, the cross section is larger at 1 TeV because of the 
't-channel' vector boson propagators. 

Figure [5] shows the rate for the W ± Z°H T associated production process ( |2.5| ). 
The cross section is small for all parameter values that we are considering. This is 
because the amplitude is proportional to sin(/3— a) cos(/9— a) and this is suppressed in 
the MSSM as, in the decoupling limit when M H ± becomes large, cos(/3 — a) becomes 
small. There is also a cancellation between the h° and H° mediated diagrams. 

The situation with process ( |2.5| ) is in good contrast with (h° / H° / A°)W ± 
associated production, (^6|) - (|2.8| ), shown in figures §. Process ( |2.6|) is enhanced 
when the decay H ± — > h°W ± has a large branching ratio. This occurs at low tan/3 
below the H + — > tb kinematic threshold. The cross section falls rapidly as the H + H~ 
pair production becomes unavailable, since there is coupling suppression as in process 

Processes ( |2.7| ) and ( |2.8| ) are not kinematically enhanced, but they are not cou- 
pling suppressed and are large compared to process (|2.5|). At large M H ± the neutral 
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Higgs boson masses also become large and these two channels are kinematically sup- 
pressed. 

The same situation as process (|2.5|) is found in processes ( |2.9| ) - ( |2.11|) , shown 



in figures |7[ The cross sections are small because of the coupling suppression and 
because of the cancellation between the h° and H° diagrams. 

It is interesting to compare the cross sections of the vector fusion processes ( |2.12j ) 
- fl2.14|) , shown in figures |j, against the s-channel processes (|2.6|) - ( |2.8|) . Although 



the vector fusion processes are at higher order in «ew; they are not suppressed 
compared to the s-channel processes especially when the centre-of-mass energy -s/s is 
large. The general behaviour of the cross sections are similar, except the h associated 
production processes ( [2. 12 ) which, compared to Q2.6|) , does not have the enhancement 



due to the resonant decay H ± — > h°W ± . Compared to the Z° associated production 
process ( |2.11| ), the cross sections are typically two orders of magnitude higher, when 



the neutral Higgs masses are small, partly because of the enhancement coming from 
the collinear photon which contributes to the neutral Higgs associated production 
modes, and partly, again, because of the cancellation between the h° and H° diagrams 
found in process (|2.11|) . The cross sections are small at a 500 GeV collider but the 



situation improves at a 1 TeV machine, where there are significant regions in M#± 
and tan (3 where the signal exceeds 10~ 5 GeV. 

4. Discussion of phenomenology 



We note that processes (|2.2j), ( |2.6| ), ( |2.7| ) and ( |2.8| ), as well as the Z — > bb decay 



mode of process ( [2.5p , which has branching ratio of 15.13% PI , an typically lead to 
the final state bbH ± W T , and the different resonance structures imply that there is 



little interference between these processes. This indicates that, as proposed in [[11 
e~e + — > bbH^W 11 is an excellent alternative mode for charged Higgs study at LC, 
even in cases where the charged Higgs mass is more than half the collider energy. 

The exact procedure for the tagging of the final state bbH ± W z ^ requires back- 
ground simulations and is beyond the scope of this study. Here we only outline the 
procedure. 

In regions where the dominant decay mode of the charged Higgs bosons is t~v t , 
we can look at the hadronic decay of the W. If we assume that the initial state 
radiation is negligible, we can estimate the four-momentum of the neutrino. The 
exact reconstruction is not possible as r~ also carries some missing momentum, and 
it is dependent on the r~ decay spectrum. 

Although there is large background coming from the decay of the top pair at 
0(0.1 pb), the reconstruction of the mass of r~v T and other kinematic cuts to sup- 
press events that can come from top pair production should reduce this background 
to a negligible level. If the initial state radiation is not negligible the signal selection 
becomes more involved, but we note that the t~v t decay mode is relevant mainly 
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when the charged Higgs mass is below the top quark mass, and in this region a 500 
GeV collider suffices, where the initial state radiation is relatively small. 

In regions where the dominant decay mode of the charged Higgs bosons is tb (or 
bbW + ), the final state is bbbbW + W~ . We can select the leptonic mode of one of the 
W's, with the other decaying hadronically, retaining about 4/9 of the total rate. We 
then reconstruct the Ws. The analysis from then on is taxing, but the total SM 
background rate for the final state tibb is 0(a| w a|) and should be small after we 
introduce a cut to eliminate soft 6's. 



The rv T associated production mode (2.1) is interesting for our purposes only 
in the large M#± region where the pair production contribution is suppressed and 
the dominant decay mode is tb. The dominant background contribution is presum- 
ably from top quark pair production followed by the decay of one of the top quarks 
into bfv T . The background reduction procedure would rely on W and t mass re- 
construction to eliminate events with two top quarks or with two Ws, and naively 
this would reduce the background by 0(a% w ), so that the signal would be visible. 
Detailed simulations are nevertheless desirable. 

In the H ± W T associated production mode, if the decay mode is T~v T , we can 
select the hadronic decay of the W . We can again find the missing momentum and 
reconstruct the charged Higgs mass, and by filtering events where this mass is small 
we can eliminate the SM W + W~ background. If the decay mode is tb we can deal 
with the top pair background by totally reconstructing the final state and eliminating 
events with two top quarks. This is expected to reduce the top pair background by 
about «ew, hence nearly two orders of magnitude. This should enable the observation 
of the peak at the charged Higgs mass in some regions of the parameter space. We 
note that the situation in this channel is better than at the hadron collider plfl as 



the background is smaller and the final state is more clean. Finally, we note that at 
fixed centre-of-mass energy the H ± W^ i cross section peaks at M H ± ~ m t + m b [|J 
as seen above, and at fixed M H ± the cross section peaks at \fs ~ 2m t |J, which is 
the energy range that is scanned over for the SM top threshold studies. At the exact 
cross section peak, the rate is one order of magnitude greater, if it is kinematically 
accessible viz M H ± + M w ± < */s. 

In all of the above cases, the measurement of the tau [22] or the top polari- 
sation will provide confirmation of the presence of the charged Higgs. 



5. Conclusion 

We discussed the single production of charged Higgs bosons at next generation linear 
colliders, and found that the cross sections are large enough to allow us to use these 
modes to study the properties of the charged Higgs bosons. Such analyses will 
complement the usual pair production process. Above the kinematical bound for the 
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pair production process, the single production cross sections are also small, although 
the signal remains visible in some channels. 

We found that the tv t H + channel, the various channels contributing to the 
bbH ± W Zf final state, and the loop induced H ± W T channel are all useful modes for 
studying charged Higgs phenomenology. The H^W^ channel is enhanced at low 
tan j3, whereas the ru T H + channel is enhanced at large tan (3. These two are the 
most promising channels for charged Higgs study beyond the kinematic limit for pair 
production. The bbH^W^ 1 channel, which has contributions from the tbH~ process 
and from the (/i°/if /v4°)W /± if =F process, is large at both large and small values 
of tan /3, but the various channels contributing to this final state are kinematically 
suppressed and the cross section falls rapidly at large M H ±. At small M H ±, this 
mode offers an attractive alternative for studying charged Higgs boson properties, as 
there are many different channels contributing to it, which could be distinguished by 
means of kinematic cuts. The analysis of this final state can then be used as a test 
of the underlying theory. 

Our analysis has been carried out at the production level only and we have 
refrained from commenting on the possible 'discovery' or even 'detection' contours. 
The evaluation of these would require a study of the decay modes, the partially 
or fully hadronic final states, and the backgrounds. The concretisation of these 
numbers would require machine dependent detector level simulations which will serve 
to complement our analysis. 
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Appendix 

The processes ( |2.1|) and ( |2.2| ), denoted generically as: 



were calculated as follows. The charge conjugate process is equivalent upto gauge 
violating effects introduced in the treatment of the widths in the propagators. 
First, we define the helicity-dependent propagators: 
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Table 1: The coupling coefficients in processes ( |2.1| ) and (|2.2|). 



Po = Pi + P2, s = Pq is the centre-of-mass energy squared, and A^ is the Z° prop- 
agator. The coupling coefficients Q and rj are given in table [I], where we defined 
2 sin 2 6w for convenience. 
In terms of these propagators the amplitude is written: 

e 3 



x w 



Mot(A e 



xV 

sMwy/Xw x 



u 3 (m 3 tan/3P L + m 4 cot (3 P R ){}( 3 + m A ) ^ Ae A 4 (A e , \)P X - 
- E Xe A 3 (X e , \)P\{^a - m 3 )(m 3 tan (3P L + m A cot(3P R ) + 
+A 5 (A e )(2p 5 - po) • E Xe (m 3 tan pP L + m A cot (3 P R )P X W (A.3) 



Here -Pl//?, = P-/+ — (1 — /+7s)/2 are the left and right chirality projection operators. 



The electron current Ex is defined as 



El = v 2 rPx e ni. 



(A.4) 



We define 
variables: 



m 4i3 (tan/3, cot /3)/M w , and the following mass dimension — 1 

(A.5) 



P(A e , A) = - /U (A) v ^[A 4 (A e , A) + A 3 (A e , A)] 
M(A e , A) = m 4/ u(A)[A 4 (A e , A) - A 4 (A e , A)] + 
+ m 3/ u(A)[A 3 (A e , A) - A 3 (A e , A)] 



(A.6) 



Po 



Q"(A C ,A) =/i(A)[A 3 (A e ,AK-A 4 (A e ,A)K-^) + A 5 (A e )(^-^)]. (A.7) 

Here k 3 = p — p 4 and k A = p — p 3 . We introduced the notation A = —A. Let 
us suppress the index A e from here on. The dependence on the incoming electron 
chirality will be implicit. The total amplitude is now written as: 

c 3 
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(A.8) 



The requirement of gauge invariance implies that the above amplitude is equivalent 
to the amplitude due to a Goldstone boson if we replace E by the momentum carried 
by the gauge boson p , and set all widths equal to zero. This condition can be stated 
in terms of the above variables as follows: 



2Qa • Po + Pav^ = 0, 



(A.9) 
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_s-2p -P3 s-2p -P4. 
Squaring up equation (|A.8|) and summing over final state helicities, but not 
dividing by 4 for combinations of initial state chiralities, we obtain: 
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We defined p x such that p Xe = pi and p Ae = p2- We observe that the above expression 
is symmetric under the exchange (3 4), (M A <-> M A ), (Q A <-> — Q A ), as the ampli- 
tude is. This can easily be seen by running a charge conjugation operator through 
equation ( |A.8|) . For numerical evaluation we can take: 



Pi 

P2 



1,0,0,1) 



^fl, 



0,0,-1) 



so that for any 4- vector and up to the phase of E Xe we have: 

A.p x = y^(A°- \\ e A 2 
A-E Xe = y^(A x + i\ e A y ) 



(A.12) 
(A.13) 

(A.14) 
(A.15) 
(A.16) 



The evaluation is relatively straightforward, and easily testable using the gauge in- 
variance and charge conjugation symmetry tests given above. 
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(a) processes ( J2.1|j and 
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(b) process Q2.3|) . Shaded circles represent one- loop contributions. 




- -H 



(c) process 



. — if ± mixing is taken into account by using the renormalised 
couplings. We neglect the box diagram. 
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(g) processes (gl^) - (gl3D 



Figure 1: Feynman diagrams for the single production processes. In all diagrams H 
stands for and W stands for , the charge being dictated by charge conservation in 
the diagram concerned. Charge conjugated diagrams have been omitted from (b), (c) and 
(f) for simplicity. 
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Figure 2: Total cross sections for the tau and the heavy quark associated production 
channels. 
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Figure 3: Total cross sections for the W H + associated production process (|2.3| 
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(a) Process fl2.6| ) 
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(b) Process Q 
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Figure 6: Total cross sections for the neutral Higgs associated production channels. 
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Figure 7: Total cross sections for the vector fusion mediated gauge boson associated 
production channels. 
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Figure 8: Total cross sections for the vector fusion mediated neutral Higgs associated 
production channels. 
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Figure 9: Total cross sections for the pair production channel (1.1) 
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